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Abstract 

We study the vacuum state of spin chains where each site carry an arbitrary representation. 
We prove that the string hypothesis, usually used to solve the Bethe ansatz equations, is valid 
for representations characterized by rectangular Young tableaux. In these cases, we obtain the 
density of the center of the strings for the vacuum. We work out different examples and, in 
particular, the spin chains with periodic array of impurities. 
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1 Introduction 



Since the seminal work of H. Bethe [1] solving the spin 



I Heisenberg spin chain [2], numerous 



applications of this method have been appeared to solve more involved models: anisotropic XXZ spin 
chain [3-5]; spin 1 chains [6,7]; alternating spin chains [8-11]; spin chains with higher spins [5,12-19]; 
spin chains with boundaries [20-23] or with impurities [24-26]. In the papers [15,27], the Bethe ansatz 
equations (BAE) have been obtained for spin chains based on gl{N) and where each site may carry 
a different representation. This generic framework encompasses the models with higher spin on each 
site or with impurities and the alternating spin chains. 

In this paper, we are interested in the resolution of these equations within the thermodynamical 
limit, but keeping arbitrary the representation on each site. The first step consists in using the string 
hypothesis. We show that the consistency of the BAE implies a restriction on the representations 
involved in the spin chain. This restriction is solved by taking representations characterized by 
rectangular Young tableaux. For the other cases, the validity of the string hypothesis remains an 
open question. Then, we determine the vacuum state and prove that it is non degenerate and a trivial 
representation of the sl{N) symmetry algebra. This vacuum state is constructed on strings whose 
lengths and multiplicities are explicitly given. Finally, the densities of the center of these strings are 
computed exactly. Throughout the paper, different examples are worked out. In particular, we treat 
in detail the spin chains with periodic array of impurities for which we compute the energy of the 
vacuum state. 

The plan of this paper is as follows. In section [2l we recall basic notions and notations on Yangian 
algebras and spin chains. Section [3] is devoted to the study of the BAE within the string hypothesis, 
and the constraints it imposes on the representations entering the spin chain. The vacuum state is 
determined in section HI see theorem 14.11 Section [5] is concerned with the resolution of the BAE 
within the thermodynamical limit. In particular, we compute the densities of the vacuum state, see 



theorem 15. 2[ 



2 Generalities 



2.1 Yangian y{gl{N)) 

We will consider the gl{N) invariant R matrices [28,29] 



i Vab 




X 



where Vab is the permutation operator 



N 




(2.2) 



and Eij are the elementary matrices with 1 in position {i,j) and elsewhere. 
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This R matrix satisfies the Yang-Baxter equation [4, 5, 28-30] 

Rabija — K) Raci^a) Rbc{K) = Rbc{K) Rac{^a) Rabija — ■ (2.3) 

It can be interpreted physically as a scattering matrix [3,4,31] describing the interaction between 
two solitons (viewed in this framework as low level excited states in a thermodynamical limit of a 
spin chain) that carry the fundamental representation of gl{N). 

The Yangian y{gl{N)) [32] is the complex associative unital algebra with the generators T-^^^ 
{Tj^p = 6ij), for 1 < i,j < N,n G Z>o subject to the defining relations, called FRT exchange 
relation [33], 

RabiK - Xb) Ta{\a) %{\) = %{\) Ta{\a) Rab{K ' \) , (2.4) 

where the generators are gathered in the following matrix (belonging to End{C'^) y {gl{N))[[X~^]]) 

N N 

T(A) = J2 ® T,,(A) = E ® E = E ^^"^ • (2.5) 

i,j=l i,j=^ r>0 r>0 

Using the commutation relations (12.41) . it is easy to show that T^^^ generates a gl{N) algebra. 



2.2 Algebraic transfer matrix 

In the following, in order to construct spin chains, it will be necessary to deal with the tensor product 
of L copies of the Yangian. For 1 < z < L, we denote by Cai{X) G End{C^) ® y{gl{N)) one copy 
of the Yangian which acts non trivially on the z*'^ space only. The space a, always isomorphic to 
End{C^) in the present paper, is called auxiliary space whereas the space i is called quantum space. 
Mimicking the notation (12. 5p used for T(A), we will consider Cai{X) as an N x N matrix whose 
entries are operators acting on the quantum space i. These operators will be expanded as series in 
A, so that we have 

£..(A) = ® L,,{X) = E ® E 4? = E 4? • (2-6) 

Obviously, £ai(A) satisfies the defining relations of the Yangian 

Rab{Xa — Ab) Cai{Xa) ^bi{Xb) = Cbi{Xb) Cai{Xa) Rab{Xa — Afo) . (2.7) 

Let us stress that the matrix Cai{X) is local, i.e. it contains only the i*^ copy of the Yangian. On the 
contrary, one constructs a non-local algebraic object, the monodromy matrix 

Ta{X) = CaliX) Ca2{X) . . . G Enrf(C^) ® {y{gl{N))f'' . (2.8) 
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Let us remark that the quantum spaces are omitted in the LHS of (12.81) . as usual in the notation of 
the monodromy matrix. The entries of the monodromy matrix Ta{X) are given by 

TV 

T^JW= ^^'^•i(^) ® Lk.kM ® ... ® Lk,_,,iX). (2.9) 

ki,..., kL-i=l 

and satisfies the defining relations of the Yangian 

RabiXa - \b) Ta{\a) %{\,) = %{\,) Ta{\a) Rab{\a - Afe) . (2.10) 

Now, we can introduce the main object for the study of spin chains, i.e. the transfer matrix 

N 

t{\)=tra{Ta{\)) = Y,Tu{\) . (2.11) 

i=l 

Equation (12.101) immediately implies 

[t(A), t(/i)]=0 (2.12) 

which will guarantee the integrability of the models. 

Let us remark that, at that point, the monodromy and transfer matrices are algebraic objects 
(in [y {gl{N)))®^) , and, as such, play the role of generating functions for the construction of mon- 
odromy and transfer matrices as they are usually introduced in spin chain models. The latter will 
be constructed from the former using representations of the Yangian, as it will be done below. 

2.3 Representations 

In order to construct representations of y{gl{N)), the following algebra homomorphism (called eval- 
uation map, see e.g. [34]) from y{gl{N)) to U{gl{N)) (universal enveloping algebra of gl{N)) will be 
used 

L,,{\) ^ 6,, + '-^ , (2.13) 

where Cpq {1 < p,q < N) are the abstract generators of the Lie algebra gl{N). They satisfy the 
following commutation relations 

[cpg, eki] = 5gk Cpi - 5ip Ckq . (2.14) 

Spin chain models will be obtained through the evaluation of the algebraic monodromy and transfer 
matrices in Yangian representations. We thus present here some basic results on the classification of 
finite-dimensional irreducible representations of y{gl{N)). 
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2.3.1 Evaluation representations 

Keeping in mind the forthcoming spin chains interpretation, we choose for each local y{gl{N)) 

algebra an irreducible finite-dimensional evaluation representation. 

We start with a finite-dimensional irreducible representation of gl{N), M(-n7), with highest weight 
= {zu^^\ . . . , zu^^^) and associated to the highest weight vector v. This highest weight vector obeys 

ekjV = , l<k<j<N (2.15) 
Cfcfc V = w^^'^ V , l<k<N, (2.16) 

where w'^^\ . . . , w^^^ are integers with w^^^'^'' < w^^'^ (these constraints on the parameters w^^'^ are 
criteria so that the representation be finite-dimensional and irreductible) . 

The evaluation representation M\{'uj) of y{gl{N)) is built from M{'uj) and follows from the 
homomorphism (12.131) . according to 

Ljk{\) y = , l<k<j<N (2.17) 
Lkk{X)v=\l + ^-\ V , l<A;<iV. (2.18) 



The representation Ma((1, 0, . . . , 0)), associated to the gl{N) fundamental representation, of C-{\) 
provides the R matrix 02. II) . 

2.3.2 Representations of the monodromy matrix 

The evaluation representations of C{X) allow us to build a representation of the monodromy matrix. 
Indeed, evaluating each of the local Ca^i{X) in a representation Mx+ie^'^e) for 1 < £ < L, the tensor 
product built on 

M^+ieA^i) ® ■ ■ ■ ® Mx+^e^{^L) (2.19) 

provides, via 02.91) . a finite-dimensional representation for T{X). The shifts iO^ entering in the defini- 
tion of the representations are called inhomogeneity parameters. Denoting by vi the highest weight 
vector associated to tu^ = {wf^^ . . . , w^^^)^ the vector 

v+ = vi®---®vl (2.20) 

is the highest weight vector of the representation 02.191) i.ejf] 

Tjk{\) t;+ = , l<A;<j<A^ (2.21) 

L 

Tjj{\) v+ = ^(^\ + i9t + i w^p'^ v+ , l<j<N. 



[2.22] 



^We have renormalized the generators in such a way that the represented monodromy matrix is analytical in A. 
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In the following, we use the following notation, for 1 < j < 

L 

Pj{\) = \[[\ + ie, + iw'^P'^ . (2.23) 
1=1 

These polynomials, related to Drinfel'd polynomials, are usually introduced to classify the represen- 
tations of Yangians. 

2.4 Pseudo-vacuum and dressing functions 

We now compute the eigenvalues of the transfer matrix t(A). As a by-product, they will provide the 
Hamiltonian eigenvalues. 

The highest weight vector fl2.20p is obviously an eigenvector of the transfer matrix. Indeed, one 
gets 

N 

t{X) v+ = J2 Tkk{X) v+ = A°(A) v+ (2.24) 
fc=i 

with 

N 

A°(A) = P.(A) . (2.25) 

k=l 

Note that A°(A) is analytical. In the context of the spin chains, the highest weight vector v + is called 
the pseudo- vacuum. The next step consists in the ansatz itself which provides all the eigenvalues of 
t{u) from AO (A). 

We make the following assumption for the structure of all the eigenvalues of t{u) 

N 

A(A) = DkW , (2.26) 

k=l 

where Dk{X), the so-called dressing functions, have been determined in [27]. They are rational 
functions of the form 

mC'-i) , _ \{k~l) i (fc+1) Aft*) , _ \ {k) i (fc-2) 
n=l ^ "T 2 n=l ^ ^" 2 

where M(°) = M^^) = 0. The parameters M^^^ G Z+, A; = 1, . . . , A^ and aI''^ G C, n = 1, . . . , M^^\ 
k = 1, . . . ,N are determined through the celebrated Bethe Ansatz equations (see below). One can 
also relate the parameters M^''\ k = 1, . . . , N, to the Tj^^ eigenvalues, according to 

L 



T^f w = (^M(^-i) - M(^') + ^ ^£ + w?) w (2.28) 
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where w is the transfer matrix eigenvector with eigenvalue (I2.26p . 

When considering the sl{N) generators Tjj^ — tJ^\ the eigenvalues read 

L 

Sj = m(^'-i) + M(^'+i) - 2 M(^') + (2-29) 

i=i 

where a^p = - G Z>o. 

Remark that for s/(2), 5*^ corresponds to twice the 'usual' spin. 



3 String hypothesis 

3.1 BAE for the center of strings 

The following set of Bethe ansatz equations, for j = 1, . . . , N — 1 and k = 1, . . . , M^^^ has been 
obtained in [15,35] (see also [27]): 



(j) ijN a/0) mO+1) 

(i) ij • 



p (\^^> — 11 

-'i+il^fc 2 J 1=1 t=i i=i 

where we defined 



2 

Aj) _ (j) (j+i) 



= 3^ . (3.2) 



Knowing that af = w/ — wl G Z+, the L.H.S. of (13.11) can be written as 

We postulate the string hypothesis which states that all the roots {A^"'^ . . . , ^^l^i^j)} are gathered 
into Vm strings of length 2m + 1 (m G |Z+) of the following form 

Amka = ASfc + ^"' a = -m,-m + l,...,m (3.4) 

where /c = 1, . . . , Vm and a£^^, the center of the string, is real. 

The string hypothesis is sustained by numerical calculations, done for a very large number of 
types of strings. One has to keep in mind that it is valid only in the thermodynamic limit L — > +00, 
which will be always the case that we will consider below. Under this hypothesis, we have 

(2m + 1) z/^;) = M(^) . (3.5) 

Because of the vanishing of M^^^ and M^^\ the parameters v'^ and vanish also. We can 
deduce from (13. ip a set of equations for the centers of the strings, for j = 1, . . . , — 1, m G ^Z+, 



k = 1, . . . , Indeed, considering products of BAE's coming from a string, we get: 
L ( ^(i) A 

n ^ (^s. + + (3.6) 



n 



\ 1=1 «=i £=1 / 

where (p, €) ^ (m, /c) means that when p = m, one has to discard the term corresponding to £ = k 
in the product. We have introduced 

Elr\X) = ep_2™(A)ep_2„^+2(A) ...ep+2„^(A) (3.7) 
^im,n).,, /eM(A) (e4™(A)e4;n-2(A) ...e2(A)) if m = n 

(A) = < ^ 2 ^ (3.8) 

[e2(m+n+l)(A) ( e2(m+n)(A) e2{m+n-l)(A) . . . e2|m-n|+2 ( A) j e2|m-n|(A) if m 7^ TT, 

Ei7''^''(A) = (^e2|m-n|+l(A)e2|m-„|+3(A) . . .e2(m+n)+l(A)) (3.9) 

3.2 String hypothesis and constraint on the representations 

From the string hypothesis, the A^^^ parameters are real, so that the r.h.s. of equation (13. 6p is of 
modulus one. Then, the l.h.s. of this equation must also be of modulus one. This condition implies 
some constraints on the type of representations that can be on the chain: 

Proposition 3.1 The BAE of the string centers are consistent if and only if the representations 
entering the spin chain solve the equations 

Ehf + ^j =Eh?'-^) ' o<p<L, i<,<iv (3.10) 

•where we have introduced 

Proof: Demanding the l.h.s. of (13. 6p to be of modulus 1, is equivalent to ask its conjugate to be its 
inverse, that is 

L m 

n n ^.fh-^x^ + '-a^A^ - = 1 (3.12) 

1=1 q='~m 

where for simplicity we abrievated A,|^^^ into A. These equations must be fulfilled whatever the type 
of strings, i.e. for all m. Then, they are equivalent to 

t{e^p{^ + ^l?) e_„.,(A-z7l^))=l, 

£=1 
L 



>o 



i=i 



which leads to 

^ ( U) \ / U) \ L / U) \ ( 0) \ 

n + + \A ( ^- + (-r? + = n f^i + ^ - ( ^- + - \f) 

where A± = A ± Since this equation must be satisfied for all A and n, it must be obeyed for all 
A±. Then, it is equivalent to 

Annihilating the coefficients of this polynomial in A^ gives (after some algebra) equations fl3.10p . ■ 



One has to keep in mind that the string hypothesis is valid only in the L — > oo limit. Then, the 
equations (13.101) must be solved for any value of L. These equations can be alternatively written as 

E -'P E + \) - -y) = • (3-14) 

Hence, a sufficient condition to solve these constraints is to ask ci'P'')'^'^ = for all values of I and j, 
that is 

„(i) , _ A 

e,+ ' ^ = or a['^ = 0, V£ = 1,...,L, Vj = l,...,iV-l. (3.15) 

Assuming moreover that no site carries a trivial representation, i.e. V£, there exists at least one j 
such that ^ 0, the above equation is equivalent to the requirement that for all site i this j is 
unique. We will call ji this unique j and the unique non- vanishing . 
The corresponding gl{N) and sl{N) highest weights read 

gl{N) : (wP,wP,...,wP,'CC!'l^-ai,...,wf^-ai) (3.16) 



sl{N) : (0^_^,af,0^_^ (3.17) 

ji-l N-l-je 

and the inhomogeneity parameters are related to the gl{l) eigenvalue t<7^^^ through: 

9, = ^^-4'\ V£=1,...,L. (3.18) 

Once the sl{N) representation has been chosen for a site, the value of w^^^ is still free and can be used 
to fix the inhomogeneity parameter to an arbitrary value (see examples below and in section . 
In terms of sl{N) representations, this constraint means that we will consider only representation 
associated with rectangular Young tableau, see figure [H 
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an 

Figure 1: Young tableau of sl{N) 



From the above discussion, one could wonder whether other kinds of spin chains (these which do 
not solve the equations fl3.10p ) are well-defined. In fact, these spin chain models do exist and are 
integrable: the constraint fl3.10p is just a consequence of the string hypothesis, and is not related to 
the existence of the model. 

From now on, we will restrict ourself to the case of rectangular Young tableaux, so as to use the 
string hypothesis. It should be clear that this is a restriction, and that other types of spin chains 
allowing string solutions do exist. However, it does not seem to exist general classes as the one given 
by fl3.16p . We will come back on this point in section [5Jl 

Let us also remark that the above restriction permits to consider all the representation of s/(2) 
and all the fundamental representations of sl{N). It also allows to have different representations on 
each site of the chain. 

Examples Throughout the paper, three examples will be taken to illustrate the general formulae. 
We will consider: 

1. A fundamental sl{N) spin chain, that is a spin chain where all the sites carry a sl(N) funda- 
mental representation i.e. ji = 1 = a^, As explained above, the inhomogeneity parameters 
on each site may be chosen arbitrarily thanks to the gl{l) parameters wf'K For example, we re- 
cover the case where all the inhomogeneity parameters vanish if wf^ = 1 i.e. the corresponding 
gl{N) highest weight is 



2. A sl(2) spin-s chain, which is a spin chain where all the sites carry a spin s representation of 
s/(2) i.e. ae = 2s, je = l,\/£. 

3. An alternating sl{N) spin chain, where the sites 2i + 1 carry an sl{N) representation defined 
by (ji, ai) and the sites 2i carry a representation (j2, 02)- 



Due to the constraint on the representations, the equations (13. 6p reduce to equations on phases. To 
express them, we introduce 




(3.19) 



3.3 



Phases of the BAE 




and (poW = 



(3.20) 
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which obeys 
Accordingly, we define: 
$M(A) = 



^ </2p+2a(A) = ^ V52a(A) + 6'(p > 2m)ipp_2mW , p G Z+ , m G -Z+ (3.21' 

!=-m a=\^-m\+l 
2m 

V54m+2(A) + 2 ^ V'2a(A) if m = n 



a=l 



(3.22) 



n 



$5'") (A) 



V^2m+2n+2(A) + V^2| (A) +2 Yl <^2a(A) ifm^ 

a=|m— n|+l 

' 5Z <^2a+i(A) (3.23) 

a=\m—n\ 

where we have introduced the step function 6[s >r) = lifs — r>0 and otherwise. Then, the 
BAE can be rewritten as 

L 



,(j+l) 



(3.24) 



E - A^:,-^^) + E - aS) + E ^-T^(^S. - A 



p,e ) 



1=1 



where the parameters Q^^;, are integer or half-integer, depending on the type of strings. 



We shall now use the monotony hypothesis, which states that Q^'*;, increases with the Bethe's 
root This hypothesis is also confirmed by numerical calculations. As an illustration, we solved 

numerically the BAE (13.241) for the vacuum state of an alternating gl{2) spin chain with spins s = 1/2 
and s = 3/2, and a number L = 128 of sites. We then plotted the vacuum density computed in the 
thermodynamical limit (see section [S]) and its discrete versioijfl 

^(Afe) = , ^ , . , (3.25) 

where A^ are the ordered solutions. Since the densities are computed thanks to the monotony 
hypothesis, the matching between them and their discrete analogues confirms the hypothesis, at 
least for the vacuum state (see figure Hj). 

This hypothesis allows us to get the bounds on Q^^i^ sending a£\, to ±oo. A direct calculation 
shows that 

$^^'""^(±00) = ±7r (^4 min(m,p) + 2 - 5m,p^ ; $LT^(±oo) = ^vr (2 min(m,p) + l) (3.26) 
$1"") (±cx)) = ±71 min(2m + 1, a) (3.27) 

•^For gl{2), all the non-zero densities have the same curve, so that one needs not to distinguish the strings of different 
length, see below. 
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and for j = 1, . . . , N, m & 

1 ^ 

QS±oo = ±2 (^,r? + 4m + 1 + ^ 5,- min(2m + 1, a,) - J]] min(2m + 1, 2n + 1) w^i^^ (3.28) 
where we have introduced 

w'^ = 2u^£ - ut'^ - z/^+i) (3.29) 

with the convention Wm = Wm~^^^ = 0, Vm G ^Z+. 

The reached bounds are deduced from the hmiting values Q^lj^^^ shifting them by the length of 
a string: 

QSmi„ = Q£-oo + 2m+l and Qlil,,^ = qI^,]+^ - {2m + 1) (3.30) 

Since we assumed the monotony hypothesis, all the for j and m fixed, are different one from 

each other. Then, the above bounds, which indicate the possible values for the Q^ij^k, also show, 
through combinatorics, the maximal number of possible states. It is thus natural to introduce: 

Definition 3.2 For a given state, the valence of the length 2m + 1 string in the j sea is defined by 
pU) _ /nO) _ /nO) _|_ 1 fU„f „•„. 

L 

Pm - Yl { min(2m + 1, 2n + 1) w^^)} + ^ 5^^, min(2m + 1, a^) (3.31) 

where the representation at site I is given by the sl{N) weight (0, . . . , 0, a/, 0, . . . , 0) and the coefficients 

Wn^ are given in equation \3. 29\) . 

We can reformulate the spin fl2.29p of these states as 

L 



S, = - J2 (2^ + 1) ^i'^ + E ^^'n (3-32) 



We illustrate this definition with the previous three examples. 
Examples 

1. For instance, the valences corresponding to a fundamental sl{N) spin chain are given by 

P^n^ = ~ 5Z { min(2m + 1, 2n + 1) } + L (3.33) 

pO) = _ ^ |min(2m + l,2n + l)tf;(/)| for 1 < j < iV (3.34) 
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2. The valences for a sl(2) spin-s chain are 

Pm = — min(2m + 1, 2n + 1) + L min(2m + 1, s) (3.35) 

3. The valences for an alternating sl{N) spin chain are given by 

pOi) = - ^ I min(2m +l,2n+ 1) w^^'^ | + ^ min(2m + 1, ai) (3.36) 

P^J2) = - ^ I min(2m + 1, 2n + 1) w^^'^^ } + ^ min(2m + 1, as) (3.37) 

^if^ = i^^^- 5Z |min(2m+l,2n+l)u;(/)} for jVJi,J2 (3.38) 

where we have supposed ji ^ j2- 

When ji = j2 = Jo (and ai 7^ as), the valences read 

P^o) = i/^o) - ^ I min(2m + 1, 2n + 1) w^^'"^ } + ^ ( ™ii^(2"^ + 1, «i) + min(2m + 1, a2)^ 
P^J} = ^^1^- {min(2m + l,2n + !)!/;(/■)} for j V J 

4 Vacuum states 

In this section, we look for states that correspond to zero excitation (no hole) configuration, i.e. 

P^J} = ^^^\'^3 = l,---,N -I and m G (4.1) 
We call them vacuum states. We have the following theorem: 

Theorem 4.1 We consider a spin chain based on gl{N), with on each site I, a representation given 
by the Young tableau fig. [21 characterized by {a£,je). 

(i) For such a spin chain, the vacuum state is non- degenerate, and built on Vm strings of length 
(2m + 1) in the j sea, with 

L 

N'y^£ = Yl Wi,a. min(j, ji) (iV - max(j, jO) (4.2) 
1=1 

(a) The vacuum state is antiferromagnetic, i.e. it is a spin state under the sl{N) symmetry algebra 
of the transfer matrix. 
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Proof: The vacuum states are defined by Pm^ = Um , Vj = 1, . . . , — 1 and m G hence they 
must solve the equations 

L 

I min(2m + 1, 2n + 1) = ^ 6jj^ mm{ai, 2m + 1) (4.3) 

Wj = 1, . . . , N — 1 and m G 

Performing the difference of two equations (14.31) . for m and m + |, and using the identities 

min(ai, 2m + 1) = ai 9{2m > ai) + (2m + 1) 61(0; > 2m) (4.4) 

^(2m +l>ai)- e{2m > ai) = 62m.+i,a, (4.5) 

^(2m < ai) - 9{2m + 1 < a/) = 52m+i,a, (4.6) 

one gets 

L 

n>m 1=1 

This proves that the system (14. 3 p is triangular, and thus admits at most one solution, so that one 
has only to prove the existence of such a state. 
We take the values: 



m 



YSjj^5at,2m+l (4.8) 
1=1 

Pluging these values into (14. 3p . a direct calculation shows that these equations are all satisfied. Hence, 
(14.81) defines a vacuum state. 

Rewriting the expression (14.80 in terms of Um using (13.290 . we get the result (14. 2p . 

Using the values (14. 8 p and the expression (13.320 . it is then easy to compute that the sl{N) 
eigenvalues of these states vanish. ■ 



Remark 4.1 Note that the relation (14.20 implies that ^2m.+i,a, min(j,ji) (A^ — max(j, j^)) must 
be a multiple of A^. This constraint (on the existence of the vacuum state) can be viewed as a 
constraint on the length of the spin chain, with parameters a/ and ji, see examples below. 

In words, the above theorem states that the vacuum state is built with strings of length in 
each sea. The number of these strings in the sea j is a function of j and je, i = 1, . . . , L. 

Even though the construction of Hamiltonian for a general choice of representation is difficult and 
not in the scope of this work, we may expect that for an approriate choice of the coupling constants 
the spin chain becomes a bipartite lattice {Lq > 2). In this case, the Marshall's theorem [36,37] 
can be applied and the state found in the theorem 14.11 is the ground state for the corresponding 
Hamiltonian. 
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Examples 

1. For instance, the vacuum state corresponding to a fundamental sl{N) spin chain is given by 
^U)^L{N_-j)_ z.0) = o, VnG^Z>o, 1 < J < iV - 1 . (4.9) 



One recovers the usual antiferromagnetic ground state, with only real Bethe roots. Their total 

N-j) 
N ' 



number is '^'•^ ■'^ , so that the state exists only if L{N — j) can be divided by A^, i.e. L multiple 



of N. 

2. For a sl{2) spin-s chain we get 

= — and Un = , Vn G -Z , n ^ s — - . (4-10) 
We get \ strings of length 2s, as expected for the ground state. One also recovers that L must 



L 

be even. 



3. If one considers an alternating sl{N) spin chain, one gets for n G |Z 

fli — 1 a2 — 1 



= when n ^ ni = and n ^ n2 = . 



If one supposes furthermore that Oi 7^ 02, 



for i < jk 



^S={ ..f!! k = l,2 (4.12) 

^^^^^^ for J>J, 

5 Thermo dynamical limit 

The Bethe equations cannot be solved in general however interesting features of the system can be 
obtained in the thermodynamical limit (i.e. L 00). We will need the following definition: 

Definition 5.1 A spin chain is called Lg-regular (Lq G Z+) when the types of representations en- 
tering in its definition satisfy 

je+Lo = je and a^+Lo = ae , V£ = 1, . . . L . (5.1) 

The set of integers {ai, . . . , a^} corresponds to the distinct values in the set {ae\l < i < Lq} (let us 
stress that C can be different from Lq). Finally, we introduce 

= {n^ = ^^LZl\l<a<C} . (5.2) 
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We also introduce the sets of indices defined by: 

4 = {je[l,Lo] s.t. = a„} , Va e [1,/:] (5.3) 

such that 

= E (5-4) 

e=i a=i i'eia 

The cardinal |/q,| corresponds to the multiplicity of da within a subset of Lq sites. 

Examples 

1. For a fundamental sl{N) spin chain, we recall that ag = ji = 1 and thus one has Lq = 1 = C, 
U = {0} and Ji = {!}. 

2. For a sl{2) spin-s chain, we get ai = 2s, ji = 1 and one has Lq = 1 = C, J\f = {s — and 
/i = {l}- 

3. For an alternating sl{N) spin chain, one has Lq = 2 = £ if ai 7^ 02 (whatever the two values 
ji and j2 are); while Lq = 2 and £ = 1 if Oi = 02 and ji 7^ j2- We have J\f = {^^^, and 
Ji = {!}, /2 = {2} in the first case; J\f = {^^^} and Ji = {1, 2} in the second case. 

From now on, we assume that the spin chain is Lg-regular. Then, adding the constraint on the 
existence of a vacuum state (see remark 14. ip to the above regularity condition, one is led to take 
L=pNLo, p e Z+. 



5.1 Regularity and constraint on representations 

The regularity hypothesis is also a natural assumption to solve the constraint fl3.10p within the 
L — 00 limit. Indeed, this constraint applied to a regular spin chain is equivalent to 

Lo / 0) \ 2P Lo / (j) \ 2P 

E(7i^' + ^j =E(7?'-^^j . l<P<io (5.5) 

which is not affected by the limit. 

As already remarked, in the case of gl{2), the representations f l3.16p describe all the representa- 
tions, so that there is in fact no constraint for this algebra (whatever the value of Lq). Then, it is 
natural to wonder if other kinds of representations can appear when N > 2. 

When Lq = 1, a direct calculation proves that the constraint (15.50 is equivalent to the equation 
f l3.15p . so that the representations described by fl3.16p are the only ones compatible with the string 
hypothesis (whatever the value of A^) . 
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For Lq = 2, one can solve directly the equations ( I5.5p . One finds that a/ > and must fulfil 
one of the conditions (for j = 1, . . . , — 1): 

'ajj)7p) = and 4^'^ 7^^'^ = 
or 

< aj^') = 2 e^^^ 7^^') and a^^^ = -2 e^^') 7^ , e^^') = ±1 , (5.6) 
or 

/^^ = ai'^ and 7^ = -7;^') . 

Remark that the first line (when applied for all j) corresponds to the representations we have studied 
in the present paper. Plugging the value (13.111) of 7^"''*, one gets 

'a? (2e, + w'f^ + -j^ =0 and a^^^ (262 + w^'^ + w!f^'^ - = , 

or 

< 202 = e(^-) aS^'^ - (cu?^ + ti7?+^)) + J and 2 9^ = -e^') a^^^ - + ti7?'+')) + j , (5.7) 
or 

^ a['^ = and 2 {9^ + 62) = 2 j - + + mi'^ + wj'+^)) , 

that one needs to solve in O^, a[^^ = w'^p — ^^'''^^^ and tu^ , with the conditions a[^^ > 0. This is still 
an intriguing problem, and we just give the complete classification of solutions for = 3 (obtained 
by direct calculation). 

In the case of s/(3), and discarding the case of arbitrary rectangular Young tableaux on each site, 
we get eight classes of solutions, which reduce to four taking into account the symmetry between the 
two sites. In each case, one site is represented by a general Young tableau (fo, c), while the second 
one is rectangular (a, 0) or (0,a), but with a determined by b and c. We present in figure [2] the 
Young tableaux of the different cases. As previously, the inhomogeneity parameters and the gl{l) 
parameters are constrained: 

Oi + w^P = \{b+l) + c+l and 92 + w^2^ = hb+l). (5.8) 
5.2 Limit of the Bethe equations for the vacuum 

For vacuum states, using the values fl3.30p of Q,^^jnin Qm,ma.K and the monotony of the Q's, one 
gets 

Q% = k-l{i^^£ + l) k = l,...,u^. (5.9) 
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'S 



(fe,c) 



'S 



b + c+l 



(0,c) 



'S 




Figure 2: The remaining allowed representations for Lq = 2 and = 3 



Using the regularity hypothesis and equation fl5.4l) . we simplify the BAE (13.241) to: 



{n,£)7^{m,fc) 



^=1 



(5.10) 



Vfc = l,...,j.^) , Vj = l,...,Ar-l , Vme 



Note that we can restrict ourself to the cases m G A/", keeping in mind that Um = when m ^ J\f. 
In the thermodynamical limit L — s> oo, Q^^kl^ becomes a continuous variable x, whose bounds 
are obtained from (15.91) and (14. 2p : 

- Yji^ <x< Y^^ with rif) = — — E ^2m+i,a, min(j, J,) (a^ - max(j, j,)) . (5.11) 

Then, the thermodynamical limit of the BAE reads 



Y, 



(j-i) 



Yr 



0-1) 



rfy $5'"HA(^)(a;)-A(^-^)(y)) 



Yr 



U) 



Y 



dy $r'"HA(^)(x)-A(^)(y)) 



+ r^^/y$5'"^(A^^^(x)-A(f-+^^(y))| = -27rx + -lE(E^^->)*t^(^"^^(^)) (5-12) 
Vx G [-Fj/') , F^^')] , Vj = 1,...,A^-1 , VmGAT 
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(5.13) 



Making the change of variables x G [—Ym-^ , Ym^] Xm{x) g] — oo , c)o[, we get 

^ POO POO 

E / rfAar^)(A)$Lr^(Ao-A)+ / rfA ai^)(A) ^r^^Ao - A) 

+ rrfAa^^+^)(A)$L7")(Ao-A)| = -27ra;(Ao) + f J] ( 5,,,>t)(Ao) 

VAo G ] - oo , oo[ , Vj = 1, . . . , - 1 , Vm G TV 
where Aq stands for Am^(a;), and = dx/dXn\ 



5.3 Calculation of the densities 

To solve these equations, one performs a differentiation w.r.t. Xin and a Fourier transform. One 
computes: 



vi/M(A) 



rf<l>LT^(A) 
dX 

dX 

_ d^^r\x) 

dX 



m+n 



-2 E 

s=|m— n| 



's + i)2 + A2 



m + 71 + 1 



+ 



m — n\ 



m+n 



(m + n + l)2 + A2 (m-n)2 + A2 



+ 2 E 



dm— nl+l 



S2 + A2 



E 



s=|f-m|+l 

We normahze the Fourier transform as: 



A2 



2 9{- > m) 



2 -m 



(f -m)2 + A2 



fip) = ^ / e^^^ /(A) rfA . 



(5.14) 
|(5.15) 
(5.16) 

(5.17) 



Exphcitly, one finds: 

(p) 



exp ( — IpI ( max( 



m,n) 



sinh ( \p\ ( min(m, n) + |) 



sinh 



\p\ 



-2 C0sh\4 ^^T\P) - <n 



exp ( - |p|max(-,m+ -) 



^ , sinh ( |p|(min(|, m + 1^ 



sinh -y 



(5.18) 
(5.19) 

(5.20) 



We introduce \l>(p), a (A^ — 1) x (A^ — 1) matrix of £ x £ blocks: 

/I + $2(p) *-i(p) 

§_i(p) I + $2(p) ^-i(p) 



^(P) 



^_i(p) I + ^2(|>) ^-i(p) 

$_i(p) I + $2(p) / 



(5.21) 
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where the blocks are defined by 



I ^(ni.ni) 



;(n2,n2)/ 



In the same way, the BAE's r.h.s. becomes 



v 



(5.22) 



(5.23) 



and the (Fourier transformed) unknown variables an\p) are gathered in the column vector 



S(p) 



S2(P) 



with Sj (p) 



V s^_i(p) / 

Finally, one is led to the following form of the BAE: 

27r$(p)E(p) = A{p) 

Using the explicit forms of the functions, the matrix \l'(p) reduces to 

/2cosh(^) -1 
-1 

$(p) =^(p)® (-^_i(p)) , Aip) 



(5.24) 



(5.25) 



2cosh(f 



v 








-1 2cosh(^ 







-1 



2cosh(§) / 
(5.26) 

with A{p) a [N — 1) X [N — 1) matrix, reminiscent of the sl{N) Cartan matrix. Let us remark that 
the matrix A{p) is independant of the type of representations. 

From this equation, one can deduce the densities an\p) by inverting the matrix '^{p): 

s(p) = ^ {A{pr' ® ( - ^^i{p)y') A(p) 

where the inverse of the matrix A{p) is given by [38] 



(5.27) 



sinh f{N — max(z, j))— ^ sinh f min{i,j 
{A{p)-% = ^ ^ ^ 



2 



sinhlA^M") sinhfM 



(5.28) 



The previous computations allow us to find a compact expression for the densities given in the 
following theorem: 
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Theorem 5.2 Let us consider a Lo-regular spin chain based on gl{N). For the vacuum, the only 
non vanishing densities are, for 1 < k < N — 1 and ria & N {1 < a < C) , 



(fc+if-l)/2 



sm 



fe/c q={\k-ji\+l)/2 cosh ( —A ) — COS ( -^TT-\ 

where the sets la have been introduced in ( 15. 



/2tx 



(5.29) 



Kir) 



For the particular case of gl{2) spin chains, the non-vanishing densities further simplify to: 



a 



C 



(5.30) 



2Lq coshyrA' 

where Sq, = tIq, + | denote the spins on the chain, and is the cardinal of I^- 

Proof: We deduce from relation (15.271) the foUowing expression for the vectors Sj(p) defined in (I5.24p : 



Up) 



(5.31) 



with 



Aa(p) 



a=l fe/a 



, bl , 

exp [ — —ma.x[aa-,ai^ 



sinh (^Y{'^^^{(^a,(ii, 



sinh Y 



\P\ - - 
exp ( —max{aa,ac] 



sinh (^^{miia.{aa, ac] 
sinh 



(5.32) 



Remarking that the action of the matrix — \l/„i(p) on the elementary vector G (with 1 in 
position a and elsewhere) gives the vector Aq(p), we deduce that 

(-$_i(p))"' A„(p)=e„, V« = l,.../:. 

Therefore, the projection on the elementary vector of relation (I5.3ip gives 

\P\- 



( \P\\ ( \P\ 

1 sinh {N - max(2, j^))— sinh min(i, j^) — 

^e(p) = ^E ^ 1) \ 2 



2t\U 



talc 



sinh N 



]p\ 



sinh 



bl 



(5.33) 



Computing the inverse Fourier transform with the method of residues, and using the trigonometric 
linearization formula 

(a-l)/2 

sm ax sr^ 

y cos2gx , 

<?=-(a-l)/2 

one is led to the equation fl5.29p . 

The expression f l5.30p is then straightforward when focusing on then gl{2) case. ■ 
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5.4 Examples 

In this subsection, we treat the thermodynamic hmit for different examples: 

1. For instance, equation fl5.33p corresponding to a fundamental sl{N) spin chain simplifies to: 

^^(P) = T'V", .JiIk' l<k<N-l (5.34) 



sinh ((iV- A;)M) 
2tt sinh (iV^) 



and the real roots densities aQ''\X) are given by, see (15.291) . in accordance with [35,38] 



" iV cosh (f A) - cos fc) ^^-^^^ 
This result is illustrated in figure [5] where the densities (I5.35P are plotted for = 5. 

2. For a sl{2) spin-s chain, one obtains the only non-vanishing density: 

^,„i(A) = - (5.36) 

2 cosh(7rA) 



2 



3- A. If one considers an alternating sl{N) spin chain, with ai 7^ 02, the non- vanishing densities take 
the following form (we remind that = |(aa "I))) for A; = 1, ... — 1, 

. ik+M-l)/2 sin f ^) 

a«(A) = — y ^ ^ ^ ^ ^ , a = 1,2 (5.37) 

9=(|fc-ia|+i)/2 cosh A j - cos j 

in accordance with the results of [9]. 

For the particular case of an alternating sl{2) spin chain with spin Si and S2, the two non- 
vanishing densities simplify to: 

a,^_i(A) = (7,^_i(A) = 1—— . (5.38) 

'2 '^a 2^ ^ 4cosh(7rA) 

We remark that the densities do not depend on the spin, i.e. they have the same expression, 
whatever the length of the string is. Their curve is plotted in figure HI where, as a by-product, 
the comparison with the numerical solutions confirms their independance from the length of 
the string (see also footnote [3]). This fact has been noticed in [8] for an alternating XXZ spin-|, 
spin-1 chain. 

3-B. If one considers an alternating sl{N) spin chain, with oi = a2 and ji 7^ j2, the non-vanishing 
densities read (with n = |(ai — 1) and k = 1, . . . , N — 1) 



2 (fe+if-l)/2 



1 - V— -II- sm I — — I 

A^ / 



£=1 q=(\k-ji\+i)/2 cosh yj^^j ~ "^os 



In figure [HI we plot the above densities for ji = 1 and j2 = 3 with A^ = 5. There, we did not 
specify the values ai = 02, since, as above, the expression (15.391) is independent of them. 



21 



5.5 Spin chains with periodic array of impurities 

We consider a gl{N) spin chain in fundamental representations on each site except for the sites pLo 
{1 < p < L/Lo) which carry a representation with the Young tableau characterized by (j, a). These 
sites can be interpreted as periodic impurities spread along the chain [25]. In our framework, this 
chain is a particular choice of a Lo-regular chain with 

A/'={0,^} , /i = {l,...,Lo-l} and h = {Lo} . (5.40) 

We will assume that Lq > 3, the cases Lq = 1, 2 have been already treated in the previous examples. 

In this case, a local integrable Hamiltonian can be easily found amongst the conserved quantities 
in the expansion of the transfer matrix (12. lip . We fix the gl{l) eigenvalue demanding that tJj^^^ is 
such that = 0, V£. Then, the Hamiltonian can be written explicitly as 

d L L 

i7 = — ln(t(A))[^^cx ^n-i,n+ (^P-w^P+i,p(0)+O ^;-i,p(0) (5.41) 

71 = 2 p=l 

where the site is identified with L and Ca,pW = J2j k ® (A + iSkj)- The matrices £jk are the 
generators of gl{N) in the representation (j, a). 
The energy per site is 

(5.42) 

iL A(0) ^ ' 

where A (A) is given by fl2.26p and the Drinfeld polynomials are computed using the formulas fl2.23p . 
f l3A6|) and flSTTSj) : 

Pi(A) = (\ + if~TE(\ + il±^)i; (5.43) 

(X + i^-^)^ for l<k<j 
P.(A) = A--^x<; ^.2 (5.44) 

(A + i — - — ) ^0 tor k > J 

The factor i in the definition of the energy allows us to obtain a real energy. 

For the vacuum, the energy per site can be computed using the string hypothesis 

(1) 

n=l \^0,n) ' 4 71=1 \^a^ ,„J ' 4 



2 



where C = — 1 + — i^^ijj^a) ^^"^ multiplicities are given in (14. 2p . In the thermodynamical 
limit, we get 

E-C = |_^rfA(a«(A) ^ + a« (A) (5.46) 
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where the densities are given by 



sin f ^ ) 1 sin 



Mil (-) 

(i).,x_ ^o-l ^'''ViV/ (1) (\\-^ lAlZ (^A7-] 

cosh (-A j - cos (- j cosh (-A j - cos (- j 

Following the lines of [38], we transform this relation as follows (using Plancherel's theorem) 

U 7-00 smh(A'i/2) L„/_„ smh(A'x/2) ^ ' 

Note that the integrand is even. We define the new variable y = e~^^ which allows us to write 



E-C = 2^-— dy'— + — -/ dy'— (5.49) 

NLq Jq 1-y NLo Jo 



Using the following formula (see e.g. [39]) 

'i-y 

we obtain finally 



dy ^- ^— = - ^Pi^i) , (5.50) 



^ - c 4 {m - ^(^)) + ^ {H^-^ + 1) - ^(^) - m + H^)) (5.51) 

where ^{x) = ^lnr(x) is the Euler digamma function. The first term corresponds to the case of a 
chain without impurity (a = j = 1), as computed in [38], while the second term is the correction due 
to the impurities. 

On figure [3], we represent the energy per site, E — C, for different types of impurities characterized 
by a and j (with Lq = 5 and = 5). 
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Figure 3: Energy E — C in terms of a and j (for Lq = 5, = 5) 
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Figure 4: Numerical vs analytical densities for gl{2) alternating (s = |,s = |) spin chain 



25 




Figure 5: Vacuum state densities for gl{5) spin chains with fundamental representations 
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